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摘 要 :采用 数值 模拟 的 方法 ,讨论 和 对 比 了 平 壁面 固定 波形 壁面 和 运动 波形 壁面 上 的 双 马 赫 反 
射 中 的 激 波 结构 和 流 场 特点 。 计 算 结果 表明 ,固定 和 运动 波形 壁面 对 气流 的 压缩 效应 使 其 在 壁面 
处 加 速 和 减速 时 马赫 杆 发 生 明 显 的 变形 和 卷曲 。 壁 面 的 压缩 在 波 后 会 产生 曲面 激 波 和 脱 体 激 波 ， 
使 得 上 游 流 场 中 的 激 波 相互 作用 更 为 复杂 。 与 平 壁面 相 比 ,固定 波形 壁 上 的 压力 和 密度 明显 升 高 ， 
而 运动 波形 壁 上 的 压力 和 密度 随 运动 频率 的 增加 而 减弱 。 
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A numerical study on the double Mach reflection over wavy 
and traveling-wave walls 
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Abstract: In this paper, double Mach reflection over wavy and traveling-waves is investigated via numeri- 

* cal simulation. The characteristics of the shock structures and the field are presented. Results show that the 
compression and expansion induced by the curved walls induce obvious deformation of the Mach stem as 
the flow accelerates and decelerates near the wall. The curved and detached shock waves induced by the 
curved wall interact with the upstream flow , making the flow field much more complicated. Compared with 
the flows over flat wall,the density and pressure rise after shock are clearly increased, while those over 
moving walls decrease as the oscillation frequency increases. 
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激 波 反 射 问题 作为 气体 动力 学 中 的 一 个 重要 研 ” (具体 条 件 参 见 文 献 [3]) 时 ,会 形成 双 马 赫 反 射 。 
究 内 容 , 在 工程 中 有 十 分 广泛 的 应 用 后。 根据 特点 。” 它 是 由 主 激 波 、 两 个 马赫 杆 、 两 个 反射 激 波 和 两 个 接 
的 不 同 , 激 波 反射 可 以 分 为 规则 反射 和 非 规则 反射 ， 。” 触 间断 交 于 两 个 点 组 成 的 激 波 系 。 这 一 模型 在 上 世 
后 者 又 可 以 进一步 分 为 马赫 反射 和 弱 激 波 反 射 "”。 ” 纪 已 经 经 过 了 大 量 的 理论 .实验 和 数值 研究 “ 。 近 
本 研究 中 主要 讨论 双 马 赫 反 射 。 些 年 的 研究 主要 关注 一 些 更 为 复杂 的 问题 。 

当 正 激 波 以 一 定 速度 扫 过 一 定 角 度 的 斜面 倾角 。 SCHMISSEUR 等 “对 流向 和 展 向 都 有 折 角 的 平面 上 
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的 双 马 赫 反 射 问题 进行 了 研究 ,发 现 展 向 折 角 使 得 
激 波 的 压缩 性 减弱 , 激 波 后 亚 声速 区 域 消失 。 杨 肠 
ax 研究 了 三 维 垂直 双 攀 面 上 的 双 马 赫 反 射 现象 ， 
发 现 三 维 效应 会 产生 三 波 / 四 波 等 复杂 结构 。BAI 
等 号 给 出 了 马赫 反射 中 反射 激 波 和 滑 移 线 的 解析 
形式 ,并 提供 了 可 以 简单 求解 以 上 解析 解 的 代数 方 
程 。 KRASSOVSKAYA 等 (站 讨论 了 内 四 边界 上 的 激 
波 反射 问题 ,发 现 反 射 激 波 的 演化 由 反射 杭 面 决定 : 

即使 是 相同 的 激 波 马 赫 数 ,不 同形 状 和 角度 的 棉 面 
也 会 产生 不 同 的 反射 激 波 。 针 对 双 马 赫 反 射 的 马赫 
杆 的 弯曲 问题 ,SHI ^60 采用 数值 模拟 的 方法 讨论 
了 黏 性 和 高 温 气 体 效应 的 影响 , 发现 这 两 种 因素 会 
导致 射流 强度 发 生变 化 ,进而 影响 马赫 杆 的 弯曲 程 
Br. 进一步 ,SHI 等 (外 从 理论 和 数值 模拟 的 角度 讨 

ie PS IEURUM RUE OE LB A38] DET 25 tl 5 n], 3f 
IUE MARISAE BEERA. 

本 研究 为 对 双 马 赫 反 射流 动 的 主动 和 被 动 控制 
的 个 简单 的 尝试 ,希望 通过 简单 的 变形 壁面 或 运 
动 呈 面 实现 简单 的 流动 控制 ,以 改变 激 波 后 的 压力 
束 筷 。 激 波 控制 的 方法 有 很 多 ,包括 在 壁面 上 施加 
吹 艳 的 扰动 .壁面 的 水 平 震荡 等 55 。 本 研究 采用 了 
蓉 月 运动 壁面 的 策略 ,希望 通过 对 激 波 后 流 场 施加 
MESE c atc : 流 场 特性 的 目的 。 

一 由 以 上 简要 的 文献 综述 可 知 ,目前 复杂 壁面 上 
的 双 马 赫 反 射 问题 逐渐 成 为 了 研究 热点 。 而 作为 简 
单 的 流动 控制 的 模型 59 ,波形 壁面 和 运动 壁面 上 的 
双轨 赫 反 射 仍然 是 未 解决 的 问题 。 沿 无 限 长 波形 辟 
面 的 二 维 定常 流动 问 是 题 是 阿 克 莱 经 典 问题 "1!。 理 
论 钟 算 表 明 ,在 亚 音速 情况 下 , 流 线 振荡 与 壁面 协同 
一 致 ,在 各 个 方向 上 均匀 传播 ,而 在 超 音 速 流 的 情况 
下 ,壁面 扰动 沿 马赫 线 偏向 下 游 方向 传播 。 特 别 是 
在 有 激 波 的 情况 下 ,迎风 面 会 产生 激 波 , 压 力 增加 ， 
背风 面 流 道 扩张 ,压力 减 小 ,总 体 效 果 是 产生 激 波 阻 
力 。 这 一 结论 在 SUN 45! 和 TYSON 4&7! Bois jf 
直接 数值 模拟 中 得 到 了 验证 。 

本 研究 采用 数值 模拟 的 方法 ,讨论 二 维 波 形 壁 
面 和 运动 壁面 上 的 双 马 赫 反 射 问题 。 在 单个 马赫 数 
和 和 斜面 倾角 下 , 主要 分 析 波 形 壁 面 对 激 波 后 流 场 结 
构 的 影响 ,以 及 马赫 杆 的 变形 在 不 同 相位 上 的 差异 。 
本 研究 主要 内 容 如 下 :首先 ,对 本 研究 的 数值 方法 进 
行 简要 的 介绍 ;其 次 ,对 波形 壁面 和 运动 壁面 上 的 双 
马赫 反射 问题 的 计算 结果 进行 分 析 和 讨论 ;最 后 ,对 
本 研究 进行 总 结 和 展望 。 
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1 数值 方法 


本 研究 求解 的 控制 方程 是 动 网 格 下 完全 气体 的 
守恒 型 二 维 积分 型 Euler 方程 , 即 


0 
-一 F- 三 1 
RUER ndS = 0 (1) 
其 中 
plu -u'n 
pulu -u )'n+pn, 
Q- P Fens , (2) 
pv pi(u-u):ntpn, 
E : : 
p pECGu — u,) n +pu'n 


式 中 : p 为 流体 的 密度 ; u 为 两 个 方向 的 速度 ; u, 为 
网 格 运 动 速度 ; n 为 边界 上 的 法 向 量 ; p HEJ. Æ 
部 全 气体 假定 下 ,气体 满足 下 列 状态 方程 
p = DAT (3) 
-P 
pe = 2 (4) 
式 中 : y Re 和 了 分 别 为 比 热 比 .理想 气体 常数 .内 
能 和 温度 。E 为 总 能 ,定义 为 
E=e+ Qr rv) (5) 
本 研究 采用 有 限 体 积 法 求解 控制 方程 (1) 。 对 
二 维 四 边 形 结构 型 网 格 的 控制 体 , 半 离散 有 限 体积 
法 的 控制 方程 可 以 写 为 
290, + DY Fns,=0 (6) 
式 中 , Qf S 分 别 为 控制 体 体积 和 面积 矢量 。 
本 研究 的 计算 中 ,变形 网 格 的 形式 随时 间 的 变 


化 是 预 设 的 ,网 格 速度 的 计算 参考 LAT 87 中 的 方 
法 。 根 据 空间 守恒 律 


puc = > [ofu dS" a 
(1 - 6) ffu; - as"] (7) 


式 中 : 2” 和 0Q" 分 别 为 n 时 间 步 和 n+1 时 间 步 的 


网 格 体积 ; 9 为 数值 参数 ,本 研究 中 采用 Crank- 
Nicholson 格式 , 9 = 0.5 。 定 义 网 格 面 上 的 平均 速 
Euy = [jut as" , 则 上 式 可 以 写 为 
J 
AV " : 
zu i 6u 7 +(1- 0)u;, (8) 
则 网 格 速度 由 以 下 公式 计算 得 到 
nel 一 AV, 1 n 
ee v 
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控制 单元 边界 上 的 流动 通 量 的 重 构 采用 经 典 的 
Roe 格式 ,左右 两 侧 物 理 量 的 插值 方法 采用 MUSCL 
格式 。 时 间 推 进 采 用 3 阶 TVD R-K 方法 。 计 算 程 
序 的 准确 性 已 经 在 平 壁 面 上 的 双 马赫 反射 和 阿 克 莱 
经 典 问题 中 经 过 检验 ,结果 与 理论 解 符合 。 

本 研究 中 的 计算 模型 为 双 马赫 反射 ,计算 域 的 
设置 如 图 1 所 示 。 计 算 域 的 流向 (x*) 长 度 L， = 27 ， 
法 向 (y) 高 度 L，= 1 。 左 右边 界 分 别 为 人 口 .出 口 
条 件 , 上 边界 根据 斜 激 波 运动 设置 为 激 波 前 后 状态 。 
下 壁面 x < 1/6 区 域内 为 自由 流 ,其 他 区 域 为 不 可 穿 
透 条 件 。 初 始 条 件 为 从 x = 176 起 始 的 倾角 为 60。 
的 斜 激 波 ,对 应 于 马赫 数 M，= 10 的 正 激 波 扫 过 倾 
RHAH 60° 的 斜面 ,左右 状态 分 别 为 
p, = 8.0,u, =7.1447, 2 - 4.125, p, = 116.5 ， 
W or =1.4,un = 0.0,vr = 0.0,pr =1.0 
二 在 该 参数 下 ,特征 量 分 别 为 波 前 的 压力 .声速 和 
计算 域 高 度 。 该 计算 模型 首先 由 WOODWARD 
AED 提出 ,并 在 后 人 的 研究 中 广泛 作为 数值 方法 的 
NER DI. 
co 


zm 


激 波 前 后 状态 


Q 
EIE 


不 可 穿 透 条 件 


a meh 波 后 状态 x 
S4 后 状态 t 
© A1 计算 域 示意 图 
Crig.l 
己 计 算 中 的 网 格 在 两 个 方向 均 为 均匀 分 布 。 作为 
网 驹 无 关 性 的 验证 ,选取 了 3 组 不 同 的 网 格 密度 ,分 
别 为 600 x 100( 粗 网 格 )、1200 x 200 (中 等 密度 网 
格 ) 和 1800 x300( 密 网 格 )3 组 算 例 。 壁 面 压力 的 


对 比如 图 2 所 示 。 
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Schematic diagram of the computing domain 
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X 
图 2 网 格 无 关 性 验证 
Fig.2 Grid independence validation 


结果 表明 , 粗 网 格 下 的 计算 结果 在 梯度 较 大 的 
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位 置 不 能 准确 分 辨 压力 变化 ,而 密 网 格 下 的 结果 与 
中 等 密度 网 格 的 结果 一 致 。 因 此 本 研究 中 采用 中 等 
网 格 密度 进行 数值 模拟 。 


2 结果 讨论 


2.1 平 壁面 上 的 双 马 赫 反 射 


首先 关注 平 壁面 上 的 双 马 赫 反射 结果 。 前 人 对 
此 已 经 有 了 充分 的 研究 ,本 研究 仅 对 该 算 例 进 行 简 
单 的 描述 ,以 便 后 文 的 进一步 对 比 和 讨论 。 图 3 中 
展示 了 1=2 时 刻 的 压力 和 密度 云图 ,转换 为 相对 坐 
标 系 下 的 局 部 马赫 数 和 流 线 ,以 及 流动 结构 的 简单 


gts 
示意 图 。 
1.0 pz 540 
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(a) 压力 云图 
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(b) 密度 云图 
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(d) 激 波 结构 示意 图 
到 3 平 壁面 上 的 双 马 赫 反 射 压 力 、 密 度 分 布 和 流 场 结构 


Fig.3 Wall pressure, density and flow structures in double 


Mach reflection over a flat plate 

Kd rp Bib ale 3 55 B AB EAE S ARAS, TEE 
激 波 ( 图 3 (d) 中 的 i 线 ) 扫 过 流 场 时 ,在 其 下 方 会 形 
成 反射 激流 (r 线 ) 和 马赫 杆 (m 2X) ,三 者 交 为 一 点 。 
此 外 ,对 比 压力 、 密 度 和 流 线 图 ,可 以 看 到 在 交汇 点 
上 还 存在 一 条 接触 间断 的 滑 移 线 (s 线 ) 。 在 滑 移 线 
下 方 靠近 壁面 的 区 域 存在 一 股 射流 ,在 c 线 处 形成 
环流 ,使 马赫 杆 m RREZ ho SH 等 ”的 研究 表 
明 ,在 考虑 流体 针 性 时 ,马赫 杆 的 弯曲 程度 不 如 理想 
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流体 中 明显 。 在 非 平衡 流 中 ,马赫 杆 的 弯曲 程度 明 
显 大 于 完全 气体 中 的 马赫 杆 弯曲 。 在 激 波 扫 过 的 后 
方 还 存在 二 次 反射 的 激 波 (z2 线 ) 和 曲线 的 马赫 杆 
(m2 2R) ,以 及 较 弱 的 滑 移 线 (sr 线 ) ,三 者 同样 交汇 
在 一 点 。 和 需要 提 及 的 是 ,本 人 研究 中 的 计算 结果 与 
GLAZ 等 “和 SHI 等 ”的 实验 和 数值 模拟 结果 有 所 
不 同 ,主要 体现 在 x = xo 附近 的 压力 和 密度 变化 。 
文献 中 使 用 的 计算 模型 中 , 左 侧 为 真实 的 弯曲 壁面 ， 
而 本 研究 中 为 自由 来 流 ,因此 在 文献 [4] 和 |13] 中 
xo 附近 的 密度 和 压力 变化 不 明显 ,而 本 研究 的 结 
变化 较为 剧烈 。 图 4(a) ~ (b) 中 进一步 展示 了 不 同 
时 刻 壁 面 压力 和 密度 的 变化 ,图 4(c) ~ (Cd) 中 用 激 
波 扫 过 的 长 度 %, 对 横 坐 标 进 行 了 无 量 纲 化 。 
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Ex 
[I-II 
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co (2) 壁面 压力 随时 间 变化 
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(d) 壁面 密度 随时 间 变 化 (对 横 坐 标 进行 无 量 纲 化 ) 
图 4 平 壁面 双 马 赫 反 射 ,壁面 压力 和 密度 随时 间 变 化 
Fig.4 Double Mach reflection on the flat wall, 


wall pressure and density vary with time 
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计算 结果 表明 ,不同 时刻 的 壁面 压力 和 密度 是 
自 相 似 的 ,这 与 前 人 的 理论 .实验 和 计算 结果 也 是 一 
致 的 。 


2.2 变形 壁面 上 的 双 马 赫 反 射 


壁面 粗糙 度 是 实际 工程 问题 中 需要 考虑 的 重要 
因素 。 作 为 真实 粗糙 壁面 的 简化 模型 ,本 研究 中 仅 
考虑 正弦 形式 的 波形 壁面 ,变形 高 度 为 0.1。 由 于 
变形 是 沿 流向 (x 方向 ) 变 化 的 ,在 下 文中 主要 讨论 
流 场 随时 间 的 演化 。 需 要 首先 说 明 的 是 ,在 下 文中 
对 流 场 的 描述 是 自 右 向 左 计数 的 , 而 对 于 壁面 的 摘 
述 是 自 左 向 右 的 。 

在 t = 1.0 时 刻 ( 图 5a) ,激流 扫 过 的 壁面 的 法 
向 位 置 处 于 持续 升 高 的 阶段 ,壁面 的 变形 对 激 波 扫 
过 的 区 域 的 影响 尚 不 明显 ,主要 表现 为 壁面 的 抬升 
对 流动 的 压缩 而 产生 的 曲面 激流。 第 一 交汇 点 处 的 
马赫 杆 `. 反 射 激流 和 接触 间断 的 变化 尚 不 明显 。 第 
二 交汇 点 处 的 弱 接 触 间 断 在 强 激 波 的 影响 下 更 不 明 
显 ,而 曲面 激 波 、 反射 激 波 和 二 次 反射 的 马赫 杆 交 汇 
于 一 点 。 

fg t = 1.5 时 刻 ( 图 5p) , 激 波 扫 过 的 位 置 处 于 
壁面 下 降 的 阶段 。 由 于 壁面 处 流体 的 加 速 运动 ,在 
x c 2 的 壁面 附近 的 射流 使 主 马赫 杆 的 弯曲 更 加 明 
显 ,接触 间断 的 强度 明显 减弱 。 凸 起 壁面 引起 的 压 
缩 激 波 的 影响 进一步 抬升 ,使 得 原来 的 马赫 杆 发 生 
变形 ,激流 的 交汇 点 增加 为 3 个 ,第 一 交汇 点 与 前 文 
中 的 结果 一 致 ,第 二 和 第 三 交汇 点 则 变 为 曲面 激 波 
和 反射 激 波 、 曲 面 激 波 和 二 次 马赫 杆 的 交点 。 值 得 
注意 的 是 ,第 二 交汇 点 处 会 进一步 产生 反射 激 波 。 

在 激 波 再 次 扫 过 壁面 抬升 区 域 的 上: = 2. 0 时 刻 
(图 Se) , 流 场 结构 变 得 更 为 复杂 。 在 主 激 波 附近 
x2. 8 附近 的 主 马 赫 杆 形状 进一步 发 生变 化 。 由 
于 上 升 阶段 流体 的 减速 运动 ,壁面 射流 强度 减弱 ,使 
得 马赫 杆 从 上 至 下 形成 “ 先 内 四 后 外 上 同 ” 的 结构 。 
在 每 个 节点 处 均 存 在 一 个 较 弱 的 反射 激 波 。 

此 外 ,由 于 凸 起 ,在 x = 2.3 附近 同样 会 出 现 一 
个 曲面 激 波 , 但 是 在 此 阶段 尚 不 明显 。 在 主 交 点 处 
的 接触 间断 比 上 一 时 刻 更 明显 。 第 二 交点 处 的 反射 
激 波 进一步 增长 ,第 三 交点 处 的 激 波 与 上 一 时 刻 比 
没有 明显 变化 。 
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激 波 进一步 扫 过 壁面 下 降 阶 段 的 1 = 2.5 时 刻 
的 流 场 结构 如 图 5(d) 所 示 。 主 激 波 下 的 马赫 杆 由 
于 其 左 侧 复杂 射流 的 影响 而 发 生 卷 曲 , 上 一 时 刻 的 
节点 处 的 反射 激 波 随 着 节点 的 消失 而 减弱 。 在 上 一 
时 刻 被 增强 的 接触 间断 在 这 一 时 刻 进一步 增强 。 

第 二 个 凸 起 壁面 处 进一步 产生 两 个 相交 的 曲面 
激流 ,与 第 二 交点 的 反射 激 波 相互 作用 形成 复杂 的 
流动 结构 。 后 方 的 第 三 交汇 点 附近 的 流 场 与 上 一 时 
刻 相 比 仍然 没有 明显 的 变化 。 
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图 5 波形 壁面 上 的 双 马 赫 反 射 压 力 云图 


Fig.5 Double Mach reflection on the corrugated wall, 


pressure cloud maps 

图 6 中 进一步 展示 了 不 同时 刻 的 壁面 压力 和 密 
度 的 变化 。 由 于 壁面 变形 的 影响 ,结构 的 自 相 似 性 
消失 。 

对 比 前 文中 平 壁面 上 的 双 马 赫 反射 ,由 于 变形 
壁面 引起 的 流动 膨胀 和 压缩 使 得 激 波 两 侧 的 压力 比 
和 密度 比 明 显 增加 。 

由 以 上 讨论 可 知 ,壁面 的 弯曲 变形 对 马赫 杆 和 
激 波 后 的 流 场 存在 显著 影响 。 这 与 KRASSOVSKAY 
等 1 在 曲面 和 缓 变 棉 形 面 上 的 结果 是 类 似 的 。 
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(b) 壁面 密度 随时 间 变 化 
图 6 变形 壁面 上 双 马 赫 反 射 , 壁 面 压力 
和 密度 随时 间 变 化 
Fig.6 Double Mach reflection on the deformed wall, 


wall pressure and density vary with time 
2.3 ”运动 壁面 上 的 双 马 赫 反 射 


运动 壁面 是 流动 控制 的 重要 手段 ,目前 主要 用 
于 减 阻 和 降 噪 等 简单 的 流动 控制 问题 。 本 小 节 通 过 
对 运动 壁面 双 马 赫 反 射 的 数值 模拟 ,考察 该 种 形式 
的 控制 是 否 可 以 降低 壁面 上 激 波 前 后 的 压 差 。 壁 面 
的 运动 随时 间 正 弦 变 化 , 即 
q, = à * sin(x — wo)" sin (wt) (10) 
如 图 7~8 所 示 , 为 1=1.0,1.5,2.0,2.5 时 刻 的 
压力 和 密度 分 布 。 流 场 结 构 与 上 一 节 中 激 波 扫 过 固 
定 变形 壁面 的 类 似 , 在 此 不 再 一 一 效 述 。 值 得 注意 
的 是 ,由 于 壁面 的 周期 性 运动 ,在 流 场 中 会 出 现 向 上 
传播 的 脱 体 激 波 ,使 得 上 游 流 场 中 激 波 的 相互 作用 
更 为 复杂 ,主要 表现 为 接触 间断 变 弱 和 曲面 激 波 的 
相交 和 反射 所 形成 的 复杂 激 波 系统 。 对 比 上 一 小 节 
中 的 流 场 结构 ,可 以 发 现 ,在 每 个 周期 运动 的 过 程 中 
都 会 出 现 一 个 新 的 曲面 激 波 。 直 观 地 对 比 图 7 和 图 
8 中 的 压力 云图 ,可 以 观察 到 , 随 着 周期 运动 频率 的 
增加 ,运动 壁面 产生 的 曲面 脱 体 激 波 的 数量 增加 ,而 
强度 减弱 。 
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(a) 大 1.0 时 刻 的 压力 云图 (a) 三 1.0 时 刻 的 压力 云图 
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X Xx 
(b) 三 1.0 时 刻 的 密度 分 布 (b) 三 1.0 时 刻 的 密度 分 布 
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(e) 大 2.0 时 刻 的 压力 云图 
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(旧居 2.0 时 刻 的 密度 分 布 
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图 7 运动 壁面 上 的 双 马 赫 反射 ( w = 100 ) 图 8 运动 壁面 上 的 双 马 赫 反射 ( e = 200 ) 
Fig.7 Double Mach reflection on the Fig.8 Double Mach reflection on the traveling-wave 
traveling-wave wall (œ = 100) wall ( œ = 200) 
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图 9 中 进一步 展示 了 2 个 算 例 中 不 同时 刻 壁 面 
上 的 压力 和 密度 分 布 。 
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(d) 壁面 密度 随时 间 变 化 (o=200) 
图 9 运动 壁面 上 的 双 马 赫 反 射 , 壁 面 压力 
和 密度 随时 间 变 化 


Fig.9 Double Mach reflection on themoving wall, 


wall pressure and density vary with time 
与 图 6 中 的 固定 变形 壁面 的 压力 相 比 ,运动 壁 
面 上 激 波 后 压力 和 密度 明显 降低 ,意味 着 激 波 强度 
的 减弱 。 对 比 这 两 组 算 例 ,我 们 也 可 以 发 现 , 除 x =1 
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位 置 ( 即 变形 部 分 开始 的 位 置 ) 外 ,其 余部 分 的 压力 
和 密度 的 幅 值 均 有 明显 的 降低 。 对 于 e = 200 的 算 
例 , 流 场 的 自 相 似 性 也 更 为 明显 。 虽 然 目前 本 研究 
没有 进行 进一步 的 验证 ,但 是 可 以 推断 在 壁面 运动 
频率 足够 高 时 ,可 以 再 次 得 到 自 相似 的 流 场 结构 。 


3 总 结 和 展望 


本 研究 采用 数值 模拟 的 方法 ,讨论 了 变形 和 运 
动 壁 面 上 双 马 赫 反 射 的 激 波 结构 及 其 对 压力 和 密度 
变化 的 影响 。 主 要 结论 如 下 。 

1) 由 于 变形 壁面 对 气流 的 压缩 效应 ,在 凸 起 的 
位 置 会 形成 曲面 激流 。 气 流 在 壁面 附近 的 加 速 和 减 
速 使 得 马赫 杆 发 生 明 显 的 变形 和 卷曲 。 曲 面 激 波 与 
加 速 和 减速 的 气流 相 结 合 会 使 第 一 交点 的 接触 间断 
更 为 明显 。 

2) 运 动 壁面 对 流 场 的 影响 与 变形 壁面 类 似 , 由 
于 壁面 的 周期 性 运动 ,在 流 场 中 会 出 现 向 上 传播 的 
脱 体 激 波 ,使 得 上 游 流 场 中 的 激 波 相 互 作用 更 为 复 
杂 ,主要 表现 为 接触 间断 变 弱 和 曲面 激 波 的 相交 和 
反射 所 形成 的 复杂 激 波 系统 。 

3) 与 平 壁面 相 比 ,变形 壁面 上 激 波 后 的 压力 和 
密度 明显 升 高 ,运动 壁面 上 的 激 波 后 压力 和 密度 比 
变形 壁面 低 。 随 着 频率 的 增加 ,激流 后 的 压力 和 密 
度 逐 渐 趋 近 于 平 壁 面 的 结果 。 

需要 说 明 的 是 ,本 研究 中 的 结果 并 没有 进一步 
讨论 壁面 变形 高 度 .坡度 对 流 场 的 影响 ,以 及 激 波 倾 
角 马赫 数 、 黏 性 等 方面 对 计算 结果 的 影响 。 这 些 问 
题 将 在 以 后 的 研究 中 进一步 讨论 。 
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